Density of states of a binary Lennard-Jones Glass

We calculate the density of states of a binary Lennard-Jones glass using a recently proposed Monte Carlo algorithm. Unlike traditional molecular simulation approaches, the algorithm samples distinct configurations according to self-consistent estimates of the density of states, thereby giving rise to uniform internal-energy histograms. The method is applied to simulate the equilibrium, low-temperature thermodynamic properties of a widely studied glass former consisting of a binary mixture of Lennard-Jones particles. We show how a density-of-states algorithm can be combined with particle identity swaps and configurational bias techniques to study that system. Results are presented for the energy and entropy below the mode coupling temperature.Comment: 6 pages, 3 figures, accepted by J Chem Phy

Constant Pressure Hybrid Molecular Dynamics-Monte Carlo Simulations

New hybrid Molecular Dynamics-Monte Carlo methods are proposed to increase the efficiency of constant-pressure simulations. Two variations of the isobaric Molecular Dynamics component of the algorithms are considered. In the first, we use the extended-ensemble method of Andersen [H. C. Andersen J. Chem. Phys. {\bf 72},2384 (1980)]. In the second, we arrive at a new constant-pressure Monte Carlo technique based on the reversible generalization of the weak-coupling barostat [H. J. C. Berendsen et. al J. Chem. Phys. {\bf 81}, 3684(1984)]. This latter technique turns out to be highly effective in equilibrating and maintaining a target pressure. It is superior to the extended-ensemble method, which in turn is superior to simple volume-rescaling algorithms. The efficiency of the proposed methods is demonstrated by studying two systems. The first is a simple Lennard-Jones fluid. The second is a mixture of polyethylene chains of 200 monomers.Comment: 10 pages, 4 figure

Density of States Monte Carlo Method for Simulation of Fluids

A Monte Carlo method based on a density-of-states sampling is proposed for study of arbitrary statistical mechanical ensembles in a continuum. A random walk in the two-dimensional space of particle number and energy is used to estimate the density of states of the system; this density of states is continuously updated as the random walk visits individual states. The validity and usefulness of the method are demonstrated by applying it to the simulation of a Lennard-Jones fluid. Results for its thermodynamic properties, including the vapor-liquid phase coexistence curve, are shown to be in good agreement with high-accuracy literature data.Comment: 5 pages, 3 figures, accepted by J Chem Phy

Multicanonical Parallel Tempering

We present a novel implementation of the parallel tempering Monte Carlo method in a multicanonical ensemble. Multicanonical weights are derived by a self-consistent iterative process using a Boltzmann inversion of global energy histograms. This procedure gives rise to a much broader overlap of thermodynamic-property histograms; fewer replicas are necessary in parallel tempering simulations, and the acceptance of trial swap moves can be made arbitrarily high. We demonstrate the usefulness of the method in the context of a grand-multicanonical ensemble, where we use multicanonical simulations in energy space with the addition of an unmodified chemical potential term in particle-number space. Several possible implementations are discussed, and the best choice is presented in the context of the liquid-gas phase transition of the Lennard-Jones fluid. A substantial decrease in the necessary number of replicas can be achieved through the proposed method, thereby providing a higher efficiency and the possibility of parallelization.Comment: 8 pages, 3 figure, accepted by J Chem Phy

Low-temperature anomalies of a vapor deposited glass

We investigate the low temperature properties of two-dimensional Lennard-Jones glass films, prepared in silico both by liquid cooling and by physical vapor deposition. We identify deep in the solid phase a crossover temperature $T^*$, at which slow dynamics and enhanced heterogeneity emerge. Around $T^*$, localized defects become visible, leading to vibrational anomalies as compared to standard solids. We find that on average, $T^*$ decreases in samples with lower inherent structure energy, suggesting that such anomalies will be suppressed in ultra-stable glass films, prepared both by very slow liquid cooling and vapor deposition.Comment: 10 pages including appendices, 8 figures. Version accepted for Physical Review Material